Overcoming blockade in producing doubly-excited dimers by a single intense pulse 

and their decay 
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Excitation of two identical species in a cluster by the absorption of two photons of the same 
energy is strongly suppressed since the excitation of one subunit blocks the excitation of the other 
one due to the binding Coulomb interaction. Here, we propose a very efficient way to overcome 
this blockade in producing doubly-excited homoatomic clusters by a single intense laser pulse. For 
Ne2 it is explicitly demonstrated that the optimal carrier frequency of the pulse is given by half 
of the energy of the target state, which allows one to doubly excite more than half of the dimers 
at moderate field intensities. These dimers then undergo ultrafast interatomic decay bringing one 
Ne to its ground state and ionizing the other one. The reported ab initio electron spectra present 
reliable predictions for future experiments by strong laser pulses. 

PACS numbers: 33.20.Xx, 32.80.Hd, 41.60.Cr, 82.50.Kx 



Blockade is a general term and encompasses many phe- 
nomena in different fields investigated for a long time. 
Thereby a single particle prevents the flow or excita- 
tion of other particles [l[ . The first studies [H, Q were 
on single-electron tunneling induced by blockade inter- 
actions and date back to the late 1960's (see also 0- 
6]). Later on, blockade phenomena have been inten- 
sively studied in different realizations with atomic Ry- 
dberg gases, where the long-range inter-atomic dipolc- 
dipole interactions prevent the excitation of two identi- 
cal atoms by the absorption of two photons of the same 
energy [7j. The excitation blockade was shown to lead 
to, e.g., spectral line broadening enhancement of 

Penning ionization and non-Poissonian excitation 
probability distributions 12], |l3[ . Recently, blockade in 
ultracold homoatomic ensembles has attracted consider- 
able attention due to its potential application for quan- 
tum information processing 0, 0, BUI- It has been 
shown that the excitation of a single atom in the ensem- 
ble can block the excitation of other atoms (even when 
the atoms are more than 10 /im apart [l|) due to the 
dipolar shift of the atomic levels, which brings the driv- 
ing pulse out of resonance. 

The effect of excitation blockade in an ensemble of 
atoms is particularly pronounced when these atoms form 
a cluster where the excited states are also shifted by 
the binding interaction. Particularly suitable systems 
to study the phenomenon are rare-gas clusters since the 
bonding is relatively weak and the constituents preserve, 
to a large extent, their atomic character. At the same 
time, rare-gas clusters are easily amenable to experi- 
ments: they can be generated in a wide range of different 
sizes [HI from a few up to millions of atoms, al- 
lowing to study a transformation of electronic properties 
from an atom to the solid. However, even the first excited 
states of rare-gas atoms lie in the UV range and cannot 
be accessed by a single photon with conventional optical 
lasers utilized in the mentioned above experiments with 



cold Rydberg atoms. 

The advent of the n ew g eneration of light sources, like 
free electron lasers 
eration setups 



22, 




21] and high-order harmonic gen- 



allow one to produce ultrashort 
and intense coherent laser pulses of high frequencies. Ex- 
posed to very strong high-frequency pulses, clusters may 
absorb a large number of photons creating differently 
charged ions 24) even if the energy of a single photon 
is not sufficient to directly ionize them. At moderate 
field intensities, however, multiple excitation of the clus- 
ters dominates over its direct ionization 1251427 ■ These 



multiply-excited states then undergo an energy transfer 
in a very efficient way to produce differently charged ions. 
Such ionization mechanisms represent a particular case 
of the broader class of phenomena, known as interatomic 
Coulombic decay (ICD) [28j |. important in various fields 
ranging from physics to biochemistry 2^, 30j |. 

Obviously, the blockade mechanism will strongly influ- 
ence the efficiency of simultaneous excitation of several 
identical atoms in a cluster. Here, we investigate how to 
efficiently produce multiply-excited homoatomic clusters 
by a single intense laser pulse, and how the energy of 
several photons deposited on different species of a clus- 
ter is then transferred to its rest. In order to illustrate 
our findings, we concentrate on a show-case example of 
the double-excitation of Ne2 by two photons. On the 
one hand, Ne2 is amenable to accurate quantum calcula- 
tions, including all processes evoked by the strong pulse, 
as well as the underlying nuclear dynamics, and on the 
other hand, this example is of interest to experimentalists 
by itself. Importantly, it allows for a transparent inter- 
pretation of the underlying physics. The present results 
allow us to draw general conclusions which are valid not 
only for dimers but also for larger clusters. 

For the example chosen the process can be schemati- 
cally described as follows: 
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FIG. 1: (Color online) Presently computed ab initio PECs 
and ICD rates of the states relevant to the process Q. These 
are the Ne2 ground state, the ungerade Ne*(2p _1 3s)Ne singly- 
excited states, the gerade Ne*(2p _1 3s)Ne*(2p _1 3s) doubly- 
excited states, and the final ionic states Ne+(2p- 1 )Ne. Total 
decay rates for each doubly-excited state into all final ionic 
states are shown in the second panel from the top. 



Two photons of energy u> excite two Ne atoms in the 
dimer into their first excited states. These doubly-excited 
states are located above the first ionization threshold of 
the dimer and, thus, undergo ICD bringing one Ne to its 
ground state and ionizing the other one. Fig. Q] depicts 
ab initio potential energy curves (PECs) and total ICD 
rates for the states relevant for the process 0]) which were 
computed as described in detail in Refs. [3ll433j . At the 
equilibrium internuclear distance (i? e =3.1 A) of the Ne 
dimer, resonant population of the singly-excited states 
requires an energy of about 16.75 eV. The energy dif- 
ference between the statistically weighted average of the 
two groups of the singly- and the doubly-excited states 
computed at 3.1 A is about 16.03 eV. 

How to overcome the enormous blockade of about 
0.7 eV and populate the doubly-excited states of the 
dimer in an efficient way? One can, of course, use two 
laser pulses with different carrier frequencies of 16.75 and 
16.03 eV. Another possibility is to use very short laser 
pulse of femtosecond or even sub-femtosecond duration, 
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FIG. 2: (Color online) Final population of the doubly-excited 
states Ne*(2p _1 3s)Ne*(2p _1 3s) as function of the carrier fre- 
quency cj after the Gaussian pulses of different durations r 
and peak intensities Iq have expired. The decay of the doubly- 
excited states (ICD) has been excluded from the calculations. 
Open symbols: fixed peak intensity of 5 x 10 12 W/cm 2 and 
different pulse durations (see right legend). Curves: fixed 
pulse duration of 20 fs and different peak intensities (see left 
legend). Note that some results are shown on an enhanced 
scale as indicated by the xlO factor in the legend. 



in order to cover the required energy interval. Alterna- 
tively, a chirped laser pulse with carrier frequency chang- 
ing from 16.75 to 16.03 eV can be applied. These strate- 
gies are, however, not feasible at present experimental 
facilities. We suggest an alternative strategy of utilizing 
a single intense pulse to overcome this blockade. 

For this purpose we have performed full quantum me- 
chanical calculations on the double-excitation of the neon 
dimer by intense laser pulses. In order to compute the 
process ([I]) we combined the previously developed the- 
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oretical and computational approaches [27 
evaluate the excitation and decay processes in intense 
laser fields. The ensuing non-adiabatic nuclear dynamics 
has been calculated employing the efficient multiconfig- 
uration time-dependent Hartree (MCTDH) method and 
code 37 , 3|| . The values of the electron transition matrix 
elements for the 2p — ► 3s excitation, the 3s — > ep ion- 
ization of the excited states, and the direct two-photon 
2p — > el ionization of Ne atom were taken from the 
Refs. [IIEl. Calculations were performed for linearly 
polarized laser pulses £{t) — So g{t) cos cot with Gaussian- 
shaped envelopes g(t) = e~ l l T of different durations r, 
carrier frequencies w, and peak intensities Iq = £§/8ira. 

Fig. [2] depicts the population of the doubly-excited 
states Ne*(2p _1 3s)Ne*(2p~ 1 3s) as function of the car- 
rier frequency w after expiration of Gaussian pulses of 
different durations r and peak intensities Iq. In order 
to preserve the population in the doubly-excited states, 
the ICD transitions from them have been excluded from 
the calculations. One can see that irrespective of the 
pulse duration and the peak intensity the optimal car- 
rier frequency of the pulse which allows for a maxi- 
mal population of the doubly-excited states is equal to 
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FIG. 3: (Color online) Results of calculations on a three-level 
model system exposed to a Gaussian pulse of r = 20 fs du- 
ration and carrier frequency ui — 16.39 eV. Shown are the 
populations of the ground state GS (dotted curve), singly- 
excited state Rl (dashed curve), and doubly-excited state R2 
(solid curve) as functions of time. All parameters are taken 
from the realistic case of Ne2 considered here: Egs = 0, 
Er2 — 2uj = 32.78 eV, the dipole transition matrix ele- 
ment for the Rl— >R2 excitation is D — 0.36 a.u., while the 
GS— s-Rl excitation is taken to be twice more probable. Panel 
(a): the state Rl is resonant for the first excitation step, 
Em = w = 16.39 eV, and the system is predominantly singly- 
excited after the pulse has expired. Panel (b): the state Rl, 
Em = 16.75 eV, is off-resonant for the first and for the sec- 
ond excitation steps (detunings are 5 — ±0.36 eV). Therefore, 
Rl serves as a virtual state. The system is essentially only 
doubly-excited after the pulse has expired. The peak intensi- 
ties of the pulses are rather moderate and taken to lead to the 
same level of saturation in both panels (the remaining pop- 
ulation of the GS is about 62%) Of course, one needs larger 
peak intensity (7 x 10 12 W/cm 2 ) in the non-resonant case (b) 
compared to the resonant (1.15 x 10 11 W/cm 2 ) case (a). 



u) = 16.39 eV. This is exactly the photon energy which is 
required to access the doubly-excited states by two pho- 
tons Eiexc = 2w = 32.78 eV. Another important detail 
seen from Fig.[2]is that the optimum carrier frequency al- 
lows for a very efficient population of the target states (a 
few tens of percent) at rather moderate peak intensities 
of about 10 13 W/cm 2 , although it is detuned far from 
the resonant transition energies for the first (16.75 eV) 
and for the second (16.03 eV) excitation steps. 

This findings can be rationalized by a simple three- 
level model. Let us consider a system with a ground 
electronic level \GS) of energy 0, an intermediate state 
\R1) of energy Em, and a target state \R2) of energy 
Er%, which is populated by two photons of energy uj — 
Efajl. Fig. [3] depicts the populations of these levels by 
a Gaussian pulse of r = 20 fs as a function of time. The 
electronic properties of the model are indicated in the 
figure caption and are related to the presently studied 
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FIG. 4: (Color online) Final populations of the singly- 
cashed curve) and doubly-excited (solid curve) sates of Ne2 
after exposure to a Gaussian pulse of 20 fs duration and car- 
rier frequency uj = 16.39 eV as functions of the peak intensity 
of the pulse Io- The decay of the doubly-excited states (ICD) 
has been excluded from the calculations, but direct ionizations 
from the singly- and doubly-excited states which take place 
during the pulse are accounted for. The dotted curve shows 
the percentage of the dimers taken out from their ground 
state by the pulse. Note that the population of the singly- 
excited states is shown on an enhanced scale as indicated by 
the x 10 3 factor in the legend and, thus, is negligible at the 
considered peak intensities. The difference between solid and 
dotted curves indicates the impact of the direct ionizations 
which take away populations from the singly- and doubly- 
excited states during the pulse. 



Ne2- In Fig. |3ja), the state \Rl) is resonant for the first 
excitation step Em = uj. The intense pulse manages 
to transfer a large fraction of the population from the 
ground to the first excited state. During the pulse a small 
fraction of the singly-excited states is further promoted to 
the final state by the absorption of the second photon of 
resonant energy. After the pulse has expired the system 
is predominantly singly-excited. 

In Fig. [^b), both the first and the second excitation 
steps are off-resonant, but the whole two- photon tran- 
sition is resonant. In this case the intermediate state 
Em is populated only virtually, and its virtual popula- 
tion is promoted further to the target state Em during 
the pulse duration. Very importantly, after the pulse has 
expired the final population of the virtual intermediate 
state is negligibly small, and the system is essentially 
only doubly-excited. This is due to the energy conserva- 
tion law which holds for the whole transition Em, = 2w, 
but does not hold for the first step Em ^ w - Because 
of the nuclear dynamics and the large number of elec- 
tronic states participating, the presently studied Ne2 case 
is much more complicated, but the results shown in Fig. [2] 
confirm the conclusion drawn. In analogy, one may ex- 
pect for larger clusters that the target states of several 
excited atoms are optimally populated if the total energy 
of several photons fits to the energy of these target states. 

Let us now turn back to the process (p}. Fig. |4] de- 
picts the final populations of the singly- (dashed curve) 
and doubly-excited (solid curve) sates of Ne2 as functions 
of the peak intensity of the pulse after a Gaussian pulse 
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FIG. 5: (Color online) Computed electron spectra after ex- 
posure of Ne2 to a Gaussian pulse of 20 fs duration, carrier 
frequency cj = 16.39 eV, and peak intensity 1.5x 10 13 W/cm 2 . 
Upper panel: Total spectrum. This spectrum is identical with 
the ICD spectrum because the individual contribution to the 
spectrum from direct ionization (broken curve) is negligible. 
After the pulse has expired 31% of the dimers remain in their 
ground state and 69% are promoted to the doubly-excited 
states and undergo ICD (see also Fig. Lower panel: The 
breakup of the spectrum into the contributions of the four 
final ionic states (see Fig. [l|. 



with optimal carrier frequency u — 16.39 eV and 20 fs 
duration has expired. The decay of the doubly-excited 
states (ICD) has been excluded from the calculations, 
but direct ionizations from the singly- and doubly-excited 
states which take place during the pulse are accounted 
for. As was already discussed above, the population of 
the singly-excited states is almost negligible (note the 
xlO 3 factor), and at peak intensities below 10 13 W/cm 2 
all the dimers which were taken out from their ground 
state by the pulse (dotted curve) are promoted to the 
doubly-excited states. Interestingly, these dimers consti- 
tute a big fraction (69% at 1.5 x 10 13 W/cm 2 ). 

At peak intensities above 10 13 W/cm 2 ionization of 
the dimer starts to play a noticeable role. This is the 
ionization of the sing ly-excited states Ne*(2p- 1 3s)Ne by 
a second photon, and also the ionization of the doubly- 
excited states Ne*(2p- 1 3s)Ne*(2p" 1 3s) by a third pho- 
ton. These 'parasitic' processes become important at in- 
tensities of about 10 14 W/cm 2 , and considerably reduce 
the 'desired' population of the doubly-excited states (see 
the difference between solid and dotted curves in Fig. H]) . 
We also note that the direct two-photon ionization of 
the ground state is negligible at the considered peak in- 
tensities. The maximal population of the target doubly- 
excited states by the 20 fs pulse requires the optimal 
carrier frequency uj — 16.39 eV and peak intensity of 
I = 1.5 x 10 13 W/cm 2 . 



Let us now turn to the decay mechanism of the target 
states and to resulting electron spectra. The average de- 
cay lifetime of the Ne*(2p- 1 3s)Ne*(2p^ 1 3s) states at the 
equilibrium internuclear distance is around 200 fs (see to- 
tal decay rates in Fig. rjj . Consequently, ICD takes place 
after the 20 fs pulse has already expired (i.e., in a field- 
free regime) . The total electron spectrum of Ne2 exposed 
to the t = 20 fs pulse of carrier frequency u = 16.39 eV 
and peak intensity I — 1.5 x 10 13 W/cm 2 is depicted in 
the upper panel of Fig. [5] by a solid curve. At this field 
intensity, the ionization of the dimer proceeds entirely 
via the ICD mechanism, since the individual contribu- 
tion to the spectrum from the direct ionization (dashed 
curve shown on an enhanced scale) is negligible. After 
the pulse has expired, 31% of the dimers remain in their 
ground state and the rest is promoted to the decaying 
states. Therefore, the integral intensity of the spectrum 
amounts to 69%. 

Interestingly, the computed total electron spectrum 
in Fig. [3] exhibits well pronounced resonant structures, 
which are due to the vibrational levels of the final ionic 
states Ne + (2p _1 )Ne. One can learn about the individual 
ionization pathways from the lower panel of Fig. [5]where 
the breakup of the spectrum into the individual contri- 
butions of the four final ionic states of different symme- 
try is depicted. We also mention that these final ionic 
states are in part weakly bound or even dissociative (see 
Fig. rjj. Following Fig. [TJ the high-energy electrons pro- 
duced via the decays into the 2 S+, 2 II g , and partly into 
the 2 II U final ionic states result in the formation of bound 
singly-charged Ne^" dimers. The low-energy electrons as- 
sociated with the decay into the 2 E+ and partly into 
the 2 LI U ionic states result in the fragmentation of the 
dimer into a Ne + ion and a Ne atom. For larger clusters, 
several excited pairs are likely to be present and undergo 
ICD, and hence fragmentation or Coulomb explosion will 
rather be the rule. Here, it also helps that the efficiency 
of the interatomic decay grows with the size of a cluster 
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In conclusion, it is demonstrated that the blockade of 
double-excitation of identical atoms in a cluster by the 
absorption of two photons from the same pulse can be 
overcome efficiently even for rather long pulses and mod- 
erate peak intensities. High-level ab initio calculations of 
the double-excitation of Ne dimers demonstrate that the 
total energy of the two photons must fit to the energy of 
the target states, i.e. that the whole double-excitation 
process is resonant. We suggest that this rule also holds 
for the multi-photon multiple-excitation of larger clus- 
ters. This excitation scheme allows to prepare a large 
fraction of clusters in multiply-excited states, which then 
relax via interatomic decay. For Ne2 we report ab initio 
electron decay spectra which can directly be compared 
with experiment. We also note that in molecular clus- 
ters the character of the blockade mechanism and the 
efficient way to overcome it will be modified owing to 



5 



the internal degrees of freedom and manifold of partici- 
pating electronic states. Moreover, the subsequent decay 
process in molecular clusters is, as a rule, much faster 

m 



44j . Investigation of how the energy of several pho- 



tons is deposited on different subunits in a system and 
understanding of how this energy is then transferred to 
the rest of the system is an issue of general importance 
in the science of intense radiation, and the present study 
is a first step in this direction. 

The research leading to these results has received fund- 
ing from the ERC under the EU's FP7 (AIG No. 227597) 
and from the IMPRS at the MPIK, Heidelberg. 
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